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Abstract 

We suggest the procedure of direct derivation of the symmetry energy from the shift of neutron- 
proton chemical potentials AA = — Ap for nuclei beyond the beta-stability line. We observe 
the presence of anomalous strong (about 15%) shell oscillations at the symmetry energy coefficient 
^sym- Our results do not confirm the existence of exceptionally large values of the symmetry 
energy coefficient at mass number j4 ~ 100 which was earlier reported in Ref. ^. Using the 
fitting procedure, we have evaluated the volume, ?*sym,voh and surface, 6sym,surf) contributions to 
the symmetry energy. We have estimated the experimental value of surface-to-volume ratio as 
fs/v = l^sym,surf|/&sym,voi ~ 1-7 for the fitting interval A > 50. 

PACS numbers: 21.10.Dr, 21.60.-n, 21.60.Ev 
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1. The nuclear /3-stability line is derived by the balance of both the isotopic symmetry, 
i^sym, and the Coulomb, Ec, energies. However the extraction of -Esym and Ec from the 
nuclear binding energy is not a simple problem because of its complicated dependency on 
the mass number A in finite nuclei The standard procedure of extraction of the symmetry 
energy from a fit of mass formula to the experimental binding energies [2j] is not free from 
ambiguities and does not allow one to separate the symmetry energy into the volume, surface 
and curvature contributions directly. 

Moreover the symmetry energy -Egym is usually derived on the /3-stability line and some 
special efforts have to be applied to extend it beyond the ground state of nuclei [sl . On the 
other hand, new information about nuclear masses in a wide region of the stability valley 
can be used for straightforward derivation of the A-dependence of energies -Egym and Ec- 

In the present work, we suggest a non-standard procedure of extraction of the symmetry 
and Coulomb energies from the experimental data using the dependence of the isospin shift 
of neutron-proton chemical potentials AA(X) = A„ — Xp on the asymmetry parameter 
X = {N — Z)/ (N + Z) for nuclei beyond the beta-stability line. This procedure allows one 
to represent the results for the A-dependence of energies -Esym and Ec in a transparent way, 
which can be easily used for the extraction of the smooth volume and surface contributions 
as well as the shell structure. 

2. Considering the asymmetric nuclei with a small asymmetry parameter X = {N — 
Z)j 1 and assuming the leptodermous property, the total energy per nucleon EjA can 
be represented in the following form of A, X-expansion 

EjA ^eA = eo{A) + 6sym(A) + Ec{X)/A, (1) 

where Cq^A) includes both the bulk and the surface energies, 6sym(^) is the symmetry energy, 
Ec{X) is the total Coulomb energy 

Ec{X) = - — {1 - Xy (2) 
^ 20 Rc ^ ^ 

and Rc is the Coulomb radius of the nucleus. 

The beta-stability line X = X*{A) can be directly derived from Eq. ([1]) using the 
condition 



dE/A 



dX 
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where 



eciA) = 0.15Ae^/R 



c 



Along the beta-stabihty hne, the binding energy per particle is then given by 

E*/A = e*(A) + b*iA) X*' + EciX*)/A, 



(4) 



where the upper index "*" indicates that the corresponding quantity is determined by the 
variational conditions ([3]) taken for fixed A and X = X* on the beta-stability line. For any 
given value of A, the binding energy can be extended beyond the beta-stability line as 



E/A = E*/A + b:^JA)iX - X*y + AEciX)/A, 



(5) 



where AEc{X) = Ec{X)— Ec{X*). The symmetry energy bly^{A) contains the A- 
independent bulk term, &*yin,voi' A-dependent surface contribution, b^^^^^^^^A'^^^ , 



h* ( A) = b* A- b* A^^l'^ 

"symV / ^sym,vol ' ^sym.surf 



(6) 



In general, the surface symmetry energy 6* surf^ ^^"^ includes also the high order curvature 



correction ~ A 



Using Eq. ([5]), one can establish an important relation for the chemical potential Ag 
{q = n for a neutron and g = p for a proton) beyond the beta-stability line. Namely, for the 
fixed A, we obtain the following result from Eqs. ([1]) and (jl]) 



AA(X) 
where 



An — Ar 



dE 

m 



dE 
'dZ 



d{E/A) 



N 
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^K^{A) + e},{A)\{X-X*), (7) 



A. 



f^E^ 



(8) 
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On the beta-stability line, it follows Eq. ([7]) that AA(X)x=x* = 0, as it should be from 
the definition of the beta-stability line. We point out that for finite nuclei, the condition 
AA = on the beta-stability line is not necessary fulfilled explicitly because of the discrete 
spectrum of the single particle levels for both the neutrons and the protons near Fermi 
surface. 

3. The quantity d{E / A) / dX in Eq. ([7j) can be evaluated within the accuracy of ~ 1/A? 
using the finite differences which are based on the experimental values of the binding energy 
per nucleon B{N, Z) = -E{N, Z)/A. Namely, 



d{E/A) 



dX 



A 



J [B{N - l,Z + 1) - B{N + l,Z - I)] 



(9) 



Since the difference (Q is fallen for AZ = —AN = 2, the pairing effects do not affect the 
resulting accuracy. 

From the binding energy tables [l| we have obtained sets of values for AA(X, A = const) 
covering mass numbers from 8 to 238. Each set contains 3 to 11 points (X, AA/4). In Fig. 1, 
the typical sets of values are plotted for A = 100, 120 and 160. As seen from Fig. 1, the 
positions of symbols determined by ([9]) can be reproduced quite well by the hnear dependence 
on X. This allows one to extract values of quantities b*y^{A), e^i^A) and X* for a given 
mass number A, see Eq. ([7]). The linear functions AA(X) from the Eq. ([7]) obtained by the 
best fit to the experimental data are shown by solid lines in Fig. 1. 
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Fig. 1. The difference A A between neutron and proton chemical potentials as a function of the 
asymmetry parameter X for nuclei with mass numbers A = 100, 120 and 160. The experimental 
values are presented by the symbols. Solid lines show the result for AA/4 using Eq. d?]) with values 
6*yjn and e^(^) obtained from the least square fit. 
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In agreement with the Eq. ([7j), the slopes of straight lines in Fig. 1 allow us to derive 
the quantity b*y^{A) + e^^A). From the beta-stability condition AA(X) = one can also 
derive the asymmetry parameter X*{A). Finally, using the Eq. ([3]), we obtain the symmetry 
energy coefficient b^^^lA) and Coulomb energy parameter e^(y4). The corresponding results 
are shown in Figs. 2, 3 and 4. Note that the error bars in Figs. 2, 3 and 4 represent the 
standard error interval obtained for b*^^, e^{A) and X* from the least square fit for the 
corresponding set AX{X, A = const). These errors reflect the systematic deviation of AA(X) 
from the linear dependence. Experimental errors in the nuclear masses are about on order 
of magnitude less than that shown by the error bars. The error bars, however, are small 
enough to see the non-monotonic behavior (the shell structure) of the plots versus mass 
number. 
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Fig. 2. Asymmetry parameter X*{A) versus the mass number A. Filled circles - experimental 



data, dashed lines - calculations using Skyrme forces SkM and SLy230b, see Ref. 



Solid 



lines - X*{A) = 0.4A/(A + 200) (thick), and X*{A) = e*c{A)/ [bty^{A) + e^(^)] (see (H) with 
e*^{A) = 0.17^2/3, 6*yj^(A) = 26.5 - 25.6^-^/3 (thin). Open circles - beta-stability line obtained 
from the periodic system of elements. 
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In Fig. 2, we have plotted the value of X*{A) (solid dots) which was obtained from Eq. 
([7]) and Fig. 1. The dashed and dotted lines in Fig. 2 were derived from the extended 
Thomas-Fermi approximation for SkM and SLy230b Skyrme forces, respectively, see Ref. 
4]. The solid (thick) line in Fig. 2 was obtained using the empirical formula 

= Q-^^ . (10) 

The "experimental" curve X*{A) in Fig. 2 shows the non-monotonic (sawtooth) shape as 
a function of the mass number A. This behavior is the consequence of shell structure of 
single particle levels near Fermi surface for both the neutrons and the protons. Because 
of this shell structure, the Fermi levels for protons and neutrons can coincide by chance 
only. For light nuclei, this significantly affects the beta-stability line, deflecting it from the 
empirical one for stable nuclei (open circles in Fig. 2). In agreement with Eq. ([3]), the 
smooth behavior of X*{A) is achieved by a fit of the symmetry and the Coulomb energy 
coefficients. Thin solid line in Fig. 2 is obtained from Eq. ([3]) with Cq^A) = 0.17 A^^^ and 
bly^{A) = 26.5 - 25.6A-V3_ 

The A-dependence of the Coulomb energy coefficient e^(A) is shown in Fig. 3. Assuming 
Rc oc A^/^, one obtains that e^(y4) oc A^^^. Using this A-dependence of e^(y4), we have the 
best fit to the experimental values by using e^{A) = 0.17 A^^^. The corresponding result of 
the fit is shown in Fig. 3 as the solid line. The deviation of the Coulomb energy coefficient 
from the smooth A-dependence is mainly due to the shell oscillations of the Coulomb radius 
Rc. 

The dependency of the symmetry energy coefficient b*^^ on the mass number A obtained 
from the experimental nuclear masses using Eqs. ([7]) and (Q is shown in Fig. 4 as solid 
circles. This dependence shows the strong shell oscillations with the amplitude of about 
15%. For the purpose of comparison, one could recall that shell effects contribute about 1% 
to the nuclear mass. In this paper, we have performed the fit of experimental data for b*^^ 
to the leptodermous-like functional form of Eq. ([6]). To extract the expansion coefficients 
^sym.voi &sym,surf ^om the fit wc havc uscd the data with A > 12 for which the justified 
leptodermous expansion. Assuming b*-y^{A) given by Eq. as the basic dependency of the 
symmetry energy coefficient on the mass number, we have obtained &symvoi = 26.5 MeV and 
^sym.surf = "25.6 McV with the surface-to-volume ratio rs/v = l^sym,surf l/^sym.voi ~ 1- The 
corresponding function 6sym(^) is plotted as the sohd fine in Fig. 4. The calculated values of 
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the coefficients ^sym voi ^sym.surf ^ quite good agreement with the phenomenological 

ones derived from the Weizsacker mass formula 6|. The dashed and dotted hues in Fig. 4 
show the results obtained from the extended Thomas-Fermi approximation using Skyrme 
forces SkM and SLy230b, see Ref. 
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Fig. 3. The Coulomb energy coefficient e'^{A) versus the mass number A along the beta-stability 
line. Experimental data are shown as solid circles. Solid line presets Cq^A) = O.nA^^^ determined 
from the fit to experimental values. Dashed lines are the calculations using Skyrme forces SkM 
and SLy230b [41. 

One should pay attention to the sources of the uncertainty in the obtained values. Firstly, 
the shell effects in b*^^ considerably reduce the accuracy of the fit procedure. Secondly, the 
range of mass numbers covered by the experimental data does not allow one to determine 
^sym.voi unambiguously, since one cannot neglect the surface term in ([6]) even for the large 
masses {A ~ 240). Additional uncertainty is because of the small curvature corrections in 
the surface symmetry coefficient &syjn,surf- -Due to the above, one can obtain quite different 
symmetry energy coefficients 6sym,voi &sym,sm-f with approximately the same quality of the 
fit taking the different intervals of mass numbers A for the fitting procedure. In particular. 
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we have obtained fcsym,voi — 32.5 MeV, b*^^ ,,^^^ = —56.3 MeV and rs/v ~ 1.7 from the fit 
for A > 50. Earher, a similar feature of the extraction of the volume and the surface terms 



at the symmetry energy was noted in Ref. 



3[ by fitting the binding energies of nuclei. Some 



numerical results are summarized in the Table 1. Note that in theoretical calculations, the 
value of surface-to-volume ratio rs/v varies strongly within the interval 1.6 < rs/v ^ 2.8, 



see Refs. 
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Fig. 4. The symmetry energy coefficient ^gym as a function of mass number A. Symbols correspond 
to experimental data. Solid line shows b*yj^{A) = 26.5 — 25. 6^^^^^ obtained from the fit to 
experimental values using Eq. ([6]). Dashed and dotted lines present the calculations using Skyrme 
forces SkM and SLy230b, respectively (see [J]). 



Table 1. Symmetry energy coefficients b*^^^^^, ^symsurf ^^^^ surface-to-volume ratio rs/v ob- 
tained from the least square fit of symmetry energy ([6]) to b*^^ from Fig. 4 for different fitting 
intervals: A> 50 and ^ > 12. 





^sym.vob MeV 


^sym,surf) McV 


rs/v 


A > 50 


32.5 


-56.3 


1.73 


A>12 


26.5 


-25.6 


1.03 
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4. We have established the relation ([7j) between the nuclear symmetry energy and the 
isospin shift of chemical potential AA = A„ — Ap that is beyond the beta-stability line. 
This relation allowed us to evaluate the symmetry energy independently of the standard 
derivation from the mass formula. Moreover this evaluation performed for different mass 
number A, has been used to derive the "experimental" values of the A-dependent volume, 
surface and curvature terms for the isospin symmetry energy. 

The presence of shell structure in the single particle levels for neutrons and protons leads 
to the characteristic oscillations in A-dependencies of the beta-stability line, the symmetry 
energy and the Coulomb energy (nuclear Coulomb radius). We point out the presence 
of strong amplitude oscillations of the symmetry energy b*^^ which significantly exceed 
the corresponding shell effects in the binding energy. We do not observe the existence 
of exceptionally large values of the symmetry energy coefficient b*^^ for the mass number 
A ^ 100 which was earlier reported in Ref. 2|. 

Our approach allowed us to estimate the experimental values of the volume, &sym,voi) 
and the surface, &sym,surf' contributions to the symmetry energy. The result of the fitting 
procedure for &sym,voi ^sym.surf depends significantly on the intervals of the mass numbers 
A applied to the fitting procedure. By fitting all of the available information for A > 
12, we have obtained &sym,voi — 26.5 MeV and the surface contribution &*yjn,surf — "25.6 
MeV with the surface-to-volume ratio rs/v ~ 1 which is well below the corresponding 
theoretical results obtained within the extended Thomas-Fermi approximation (ETFA) with 
the effective Skyrme forces in the present work as well as in quantum Skyrme-Hartree-Fock 
(SHF) approach 10[. We also point out that the change of the intervals in mass numbers A 
for the fitting procedure leads to a significant difference in the surface contribution ftgym^surf 
to the symmetry energy, see rows 2 and 3 in the Table 1. A better agreement with theoretical 
results is obtained for the fitting interval with larger masses A> 50 where the leptodermous 
expansion is ,„ore justified. In this case, the surfptcvolume rat.o .s given i,y r.,. « 1.7 
which is close to the theoretical result from Ref. IK 
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